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1 Introduction 



Particle and transverse energy production in the central rapidity region of ultrarela- 
tivistic heavy ion collisions (URHIC) can be viewed as a combination of perturbative 
(hard and semihard) parton production and non-perturbative (soft) particle produc- 
tion. The key feature in the division is that one should be able to compute the hard 
and semihard parton production from perturbative QCD (pQCD), provided that the 
parton distributions of colliding nuclei are known |I|, 0, ||. On top of this, the non- 
perturbative particle production can be modelled e.g. through strings 0, ^, perhaps 
also through a decaying strong background colour field 0. 

The semihard QCD-processes are expected to become increasingly important with 
increasing cms-energies ^/s, particularly in URHIC, due to the following reasons: 
Firstly, already in pp{p) collisions the rapid rise of the total and inelastic cross sec- 
tions can be explained by copious production of semihard partons, "minijets," with 
transverse momenta Pt > Po ~ 1---2 GeV 0. In other words, the events tend to 
become more "minijetty" towards higher ^/s. This is also expected to happen in A A 
collisions, and even more so, because of the A^/'^-scaling of hard collisions. 

Secondly, the rapid rise of the structure function F2 observed at HERA at x ^ 0.01 
and Q^l GeV P] implies that minijet production at the LHC gets a further boost 0: 
the mid-rapidity minijets with ~ 2 GeV typically probe the gluon distributions at 
X ~ 2pt/v^ ~ 7 X 10~^ for nuclear collisions at \/s/A = 5.5 TeV. 

Thirdly, the time scale for producing partons and transverse energy into the central 
rapidity region by semihard collisions is short, typically Th ~ 1/po ~ 0.1 fm/c, where 
Po ~ 2 GeV is the smallest transverse momentum included in the computation. The soft 
processes are completed at later stages of the collision, typically at Tg ~ 1/Aqcd ~ 1 
fm/c. If the density of partons produced in the hard and semihard stages of the heavy 
ion collision becomes high enough - as will be the case - a saturation in the initial parton 
production can occur already in the perturbative region ^ Aqcd |10|, |11], |12|, and, 
consequently, softer particle production will be screened. The fortunate consequence of 
this is that a larger part of parton production in the central rapidities can be computed 
from perturbative QCD (pQCD) at higher energies, and the relative contribution from 
soft collisions with px ~ 2 GeV becomes smaller. Typically, the expectation is that 
at the SPS (Pb+Pb at ^/s = 17 AGeV), the soft component dominates, and at the 
LHC (Pb+Pb at ^/s = 5.5 ATeV) the semihard component is the dominant one. At 
RHIC (Au+Au at ^/s = 200 AGeV) one will be in an intermediate region where both 
components should be taken into account. 

The semihard processes have also been implemented in several event generators for 
URHIC [0. For example, HIJING [| and Parton Cascade Model (and VNI) [|l|] rely 



on the dominance of initial semihard particle production at high cms-energies. 

In what follows, I am going to focus on how to compute the initial minijet and 
transverse energy production in an AA-coUisions from pQCD, and what are the impli- 
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cations of such a computation for the initial conditions of QGP-formation at the LHC 
and RHIC. I will also discuss the uncertainties in the computation. 



2 Production of semihard q in AA 

Hadronic jets originating from high-px parton-parton scatterings have been observed 
in pp collisions from pt ~ 5 GeV up to p-r ~ 440 GeV |jl6| . Below 5 Ge V it 



becomes very difficult to distinguish the individual transverse energy clusters |T5| from 
the underlying background, so that minijets with px ~ 2 GeV cannot be observed 
directly Q. In AA collisions, where hundreds (RHIC) or thousands (LHC) of minijets 
are expected to be produced within the central rapidity unit, it becomes certainly 
impossible to distinguish them individually. However, due to their abundance, minijets 
are expected to dramatically contribute to QGP-formation and its further evolution at 
the LHC and RHIC. 

Unlike in a jet measurement (or next-to-leading order (NLO) calculation) we do 
not define any jet size for an individual (mini)jet, but rather, our goal will now be to 
compute minijet and transverse energy production in the whole central rapidity unit 
\y\ < 0.5. The key quantity is the integrated minijet cross section. 



(Tjetiy/s,po) = ^ fdp'^dyidy2 Xifi/Aixi,Q) X2fj/Aix2,Q) 

^ ■'Po ij<kl> = 



X 



X 



dt 



[s,t,u) 



dt 



s, u, t) 



1 + 5, 



(1) 



kl 



Collinear factorization is assumed, and fi/A{x,Q) are the number densities (per nu- 
cleon) of partons i in a nucleus A at a fractional momentum x and a factorization scale 
Q, chosen here as Q = Pt- From conservation of energy and momentum in simple 
2 — »• 2 kinematics with negligible initial transverse momenta it follows that 



Xi 



Pt 



e^^+e^^) and 



Pt , 



-yi 



(2) 



where 7/1,2 are the rapidities of the outgoing partons k, I and px is the magnitude of their 
transverse momentum. Note that in this case k and I are back-to-back in transverse 
momentum. The kinematic invariants on the parton level become 



s = X1X2S = 2p^(l + coshF), i = — Px(l + ^ 

where Y = yi — y2. 

^In this sense minijet is not a "jet" at all. 
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In the lowest order pQCD there are eight types of subprocesses a^^^^'-: 
99 99, mi 91 91 qiQi ^ QjQj, QiQi, 99 Qilj ^ Qilj QiQi ^ QiQi- (4) 



The cross sections can be found in the hterature [0. In cijet the contribution from 
99 ~^ 99 is clearly dominant. In the computation below, I will use but as a 
little side-remark it is perhaps worth noticing that the "single effective subprocess" 
-approximation |18| works also nicely: it is straightforward to show that 

[dais,i,u)/di\s^-^^ ^9^1 0{x), X « 1 

[da{s,i,u)/di\9<J^99 + [da{s,u,i)/di\9'i~*99 4 \ 0{l/x), X>1 



(5) 



where x = u/i. Considering similarly the process qq — > qq, one obtains an approximate 
relation between the different main subprocesses: 

9 9 

{99 99) ■■ i9q 9q) ■ {m^ m) = ^ - (^)^ (6) 



so that Eq. (|I]) becomes simply 



If da 99-^99 

f^jet(v^,Po) = 7T dp'^dyidy2XiF{xi,Q)x2F{x2,Q)—^ (7) 
/ Jp?, dt 



with Q) = xg{x, Q) + | Eg x[q{x, Q) + q{x, Q)]. 

Let us first neglect all nuclear effects in parton densities and take fi/A = fi/p = fi, 
which is a reasonable first approximation to start with. Then, from the pp-level quantity 
(jjet, we obtain the average number of produced minijets with > pq in an AA collision 
with an impact parameter b simply by multiplying cijet by the nuclear overlap function 
Taa 0, §: NAAih) = 2T^^(b)o-jet, where 

TAA{h) = J dhTA{s)TA{h - s), (8) 

where the thickness function is 

T4(s) = J dzuAiV^^T^) (9) 

with normalizations / d^hTAAO^) = ^ and / (i^sT^(s) = A, i.e. the hard cross sections 
in A A scale as ^4^, as expected in absence of any nuclear effects in jijA- For central 
collisions with Woods-Saxon nuclear densities ua we get Taa(& = 0) ~ ^^/(vri?^), so 
clearly A^^^ ~ A^l"^ . Here one should notice that this approach is based on independent 
binary parton collisions, and is valid when each semihard sub-collision consumes only 
a negligible fraction of the beam-energy, and when there are sufficiently many partons 
available for scattering. These conditions are fulfilled for minijet production at central 
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rapidities even up to the LHC-energies and A ~ 200. However, at the LHC for large y 
together with a very large A, one eventually runs into trouble with energy and baryon 



number conservation , which is an indication that in this regime one should consider 
coherence effects in parton production. This is easy to understand simply by noticing 
that the available energy and valence quark number scale as ~ A, while Naa ~ A'^^^. 
For RHIC energies, however, there does not seem to be a problem in the above approach 



for any y. [12 



We want to focus on the central rapidity unit, so we will have to do some book- 
keeping of where the partons are produced in rapidity, and define our acceptance cuts. 
Especially, we want to compute the average transverse energy carried by the partons 
in the central rapidity unit. Furthermore, we would like to keep track of contribu- 
tions from gluons, quarks and antiquarks separately. Let us therefore define the E^- 
distribution for each flavour / in our acceptance window Ay as 0, |T^ : 



dp%dyidy2 Yl ^ifii^i, Q) X2fj{x2, Q)^ 



^3 
{kl) 



+ 5, 



X 



kl 



X 



da'-^'^^i 

-J (s, t, u) 6{Et - [Sfke{yi) + Sfie{y2)]pT) + 
dt 



^^ij^ki ^ 

+ -J {s,u,t)5iET 
dt 



[Sfif^iyi) + 5fke{y2)]pi 



(10) 



with a normalization J2f I '^-^t||^ = crjet(-\/s, Po)- Above, our acceptance is defined 



through a "measurement function" e{y), which will in the following be chosen as a step 
function 

' 1, if \y\ < 0.5 
0, otherwise. 



:iii 



In Fig. 1, the integration region in the {yi, y2)-plane (with a fixed p^), 

\yi\ < log 



2pi 











, -log 




<y2< log 






-Pt 




-Pt 



(12) 



is illustrated together with the cuts corresponding to our acceptance Ay : \y\ < 0.5. 
After this, it is straightforward to compute the first ii^x-Kioment for each flavour / in 
the semihard collision, 



a^eti^/s,po){E^^)Ay = J dE^ET 



daf 
dE^ 



dp%dyidy2 Yxifi{xi,Q)x2fj{x2,Q) x 
1 



{kl) 



l + 5ki 



Sfk-;T is,t,u) + Sfi-j {s,u,t) 
dt dt 



PT^iVl) 



(13) 
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Figure 1: Kinematic cuts in rapidity. The oval regions show the rapidity phase space (|l2|) for 
fixed pt = 2 GeV and 5 GeV, and = 200 GeV. The values of e{yi) + £(2/2) corresponding to our 
acceptance region, Ay : |y| < 0.5, are indicated. 



In an AA collision with an impact parameter b, the average transverse energy 
carried by semihard partons of flavour / to the acceptance region becomes then |^ |l^ 

E^(b, v/s,po, Ay) = TAA{h)ajetiVs,Po){E^)Ay, (14) 

where TA^(b)cTjet(v^, Po) is the average number of semihard collisions (all y) and 
(-E't)aj/ is the average transverse energy carried by the flavour / at \y\ < 0.5 in each 
semihard collision. 

We can also compute the number distribution da^ /dn for each flavour in a similar 
way (in Eq.(|l^) replace by n, and in the delta-functions by 1) to obtain the 
average number of flavour / produced within our acceptance window: 

^iA(b, V^,P0, Ay) = TAA(b)(Tjct(A/s,Po)(^-'^)Ay, (15) 

where <yjct{y/s,Po){n^)Ay can be obtained directly from Eq. (plSf ) simply by removing 
the weight p^. The normalization I]/ crjet(v^,Po)('^"^)Ay^oo = 2(Tjet(v^,Po) can be 
verified from Eqs. (|^) and (ITSD . 
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3 Initial conditions at r = 0.1 fm/c 



The subprocess cross sections da^^ /di diverge as ~ when 0, which is why we 

cannot rehably extend our computation below po = 1...2 GeV. The QCD-evolution of 
the parton densities does not improve the situation enough because of two competing 
effects: towards smaller both x and Q decrease; the downwards QCD-evolution 
reduces the densities (for fixed, small x) but the simultaneous small-x increase of the 
densities has an effect in the other direction. We would like to extend our computation 
down to as low values of pt as possible but making sure we can still believe that we 
are dealing with perturbative quantities. Of course, only after computing the NLO 
minijet cross sections, the question of convergence of the perturbation series can be 
better answered. For large-px jet measurements the NLO-calculation has been done 
19(1 but for the few GeV px-range and for the rapidity acceptance we have in mind 



this has not been completed yet. 

So, how to fix the parameter pq7 In Fig. 2, I have plotted the behaviour of 
(^}et{y/s , Po) {E^) Ay and <Jjct{y/s,Po){n^)Ay for the LHC and RHIC energies as functions 
of Po- In the computation, I have assumed that fi/A ~ fi- I have used the GRV94 LO- 



distributions |]22|, which reproduce the observed small-x rise of the partons densities 
reasonably well. Note that there is no ad hoc "iiT-factor" included in the computation. 
One way to fix po is to formulate a saturation criterion |1|, |10|, |T2| for the produced 



parton system in URHIC: let us assign a transverse area vr/p^ for each parton within 
our rapidity acceptance. Most of the partons have pt ~ Po, so that the effective total 
transverse area occupied by the partons is Naa x tt/pq. The available nuclear transverse 
area can be estimated as ttR\, so that the saturation of parton production in \y\ < 0.5 
should occur when Naa x 7r/pQ>7ri?^, i.e. when 

/ Po \ ^ 

^jct(v^,Po)(^^)|i/|<o.5^ 134f ^^^ j mb, (16) 

where i?pb = 6.54 fm and TpbPb(O) = 32/mb has been used. This procedure results in 
Po ^ 2 GeV for the LHC, and po ^ 1 GeV for RHIC, as illustrated in Fig. 2. 

Saturation of parton production can also be studied in a more self-consistent - 
but still phenomenological - manner by computing a screening mass (electric, static) 



generated by the partons produced by a fixed time r < l/px ITlIl- This mass then 
screens further parton production at smaller px and later in time. The more partons 
are produced the stronger the screening becomes, and, the more damped further parton 
production is. In this way, saturation of transverse energy production seems again to 
take place at px ~ 2 GeV for LHC and px ~ 1 GeV for RHIC [Uj. 

Lower limits for po may also be studied in an eikonal approach 0, or, by convoluting 
the parton cross sections with the fragmentation functions of each parton flavour into 



pions and kaons pO[- Comparison of the fragmentation function approach with the 



measured px-distributions of charged particles in pp collisions is successful at ^/s > 200 
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GeV at large whereas it seems that allowing parton scatterings with < 1.5 

GeV would result in an overestimate of the charged particle px-distributions. Therefore, 
let me fix po = 2 GeV both for the LHC and for RHIC. Estimates for smaller pq can 
be easily obtained from Fig. 2. 




1.4 1.6 1.^ 

Po [GeV] 



2.2 



Figure 2: quantities <Tjct{\^,Po){n-')Ay and ajct{\/s,po){E^) Ay from Eqs. ( J13[|l^ ) for = 5500 
GeV (LHC) and = 200 GeV (RHIC) as functions of pa for Ay : \y\ < 0.5. The solid curves 
show the total results, the dashed ones the gluon contributions. For the figure, quark and antiquark 
contributions have been muhiphed by a factor 10 (7) for the LHC (RHIC) and they are shown with 
dot-dashed and dotted curves, respectively. The saturation limit from Eq. (p^ ) is the increasing 
dashed line in the upper panel. 

Table shows the average numbers and transverse energies of semihard partons at 
r = 0.1 fm/c with \y\ < 0.5 and Pt > 2 GeV in central Pb-Pb collisions at the LHC 
and RHIC, as given by Eqs. (|15D and ([l|). 

We make the following five interesting observations: 

1. Gluons clearly outnumber quarks and antiquarks and dominate E'x-pro duct ion 
at r ~ 0.1 fm/c both at the LHC and RHIC. The QGP is actually gluon plasma to a 
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\y\ < 0.5 


total 


9 






LHC: 


^^PbPb 

e{ 


4741 
14160 


4350 
12950 


200 
619 


191 
590 


RHIC: 


^PbPb 


121 
321 


99.6 
263 


13.0 
36.1 


8.10 
21.9 



Table 1: The initial conditions at t = 0.1 fm: N^^^^ and for 6 = 0, ./s/A = 5500 GeV (LHC) 
and 200 GeV (RHIC), po = 2 GeV and |y| < 0.5. 

first approximation. 

2. With po = 2 GeV ^ Aqcd the glue is saturated at the LHC. This indicates that 
the pQCD-domain is dominant for E'x-production. In other words, more gluons are 
produced into the system at p < po but their transverse energy is relatively small. For 
RHIC I expect the soft component to be still important. 

3. Estimating the volume of the system at rh = 0.1 fm/c and with Ay = 1 by 
Vi = TrRp^AyT}^, we can convert the results into local energy and number densities: 
^pQCD _ ^/^^ ^pQCD _ Nf^pjVi. From the numbers in Table I for the LHC, we 
see that E'T/gluon ^ 3 GeV (naturally, since po = 2 GeV). On the other hand, for an 
ideal gas of massless bosons ef^'^^/n^g^^^ ^ 2.7Teq. For an ideal gas with e^^^*^' = e^^'^^^ 
we get Teq ~ 1.1 GeV, which means that (eg/ng)^'^*-'^ ^ {eg/ugY'^'^'^^, and the glue 
seems to be initially thermalized at the LHC. |12| This result is essentially due to the 
small-x increase of the gluon distributions At RHIC the gluon thermalization will 
take place a little later. One should also keep in mind that a complete thermalization 
requires locally isotropic momentum distributions, and to study this more detailed 
modelling is needed. 

4. Since nq,nq <C n^, quarks and antiquarks are initially obviously far away from 
chemical equilibrium both at RHIC and the LHC. 

5. By assigning a baryon number l/3(— 1/3) for each quark(antiquark) produced 
at r = O.lfm/c, we observe that the initial net baryon number density at \y\ < 0.5 
becomes nQ_§ = 0.21 fm~^ for the LHC and 0.12 fm'^ for RHIC. This density will 
naturally dilute quite fast (~ 1/r) but it is interesting that for the LHC the initial 
number is actually larger than for RHIC, and even larger than the nuclear matter 
density 0.17 fm~^. This behaviour is again due to the small-x increase of the gluon 
densities. Also, for the LHC, it is worth noticing that the transit time of the colliding 
nuclei is 2Ra/'J ~ 0.005 fm/c ^ 0.1 fm/c. Finally, the initial net baryon number-to- 
entropy ratio can also be estimated [|1^], and we find that {B — B)/S ~ 1/5000 (1/2000) 
for the LHC (RHIC), so that we are far away from the conditions of the early Universe, 
where the inverse of the specific entropy is ~ 10~^. 
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4 Discussion and outlook 



I have shown how to compute initial conditions for QGP-production in URHIC at the 
LHC and RHIC from the lowest order pQCD. The analysis can obviously be sharpened 
in various ways: by computing the transverse energy production in NLO in and by 
including nuclear effects in the parton distributions in a consistent, scale-dependent 
Resummation similar to the Drell-Yan case [|^, if applicable, would 



manner 



hopefully help to reduce the uncertainty due to the parameter pq. Also higher twist 
effects and more coherent scattering should be studied and formulated in more detail 
at large y. Connection to the BFKL-physics potentially relevant to the minijets in- 
volving small-x should be understood better [|25| as well. Recently, a novel approach 
based on classical gluon fields has been developed for semihard parton production. 
Applicability of this model at the energies of RHIC and the LHC has been studied 
lately Also, understanding the pre-thermal evolution of the newly produced QGP 
is a challenging task, involving space-time dependent phenomena in gauge theories. 

To conclude, it is evident that there are several interesting open questions in pre- 
thermalization physics of URHIC, and much more precision-work is needed. Under- 
standing the primary parton production mechanisms will be essential for understanding 
the further evolution and signals of the QGP, and more global variables measured at 
nuclear collisions at the LHC and RHIC. 
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